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We present the first high-resolution, soft-X-ray spectrum of the prototypical Seyfert 2 galaxy, 



I NGC 1068. This spectrum was obtained with the XMM-Newton Reflection Grating Spectrometer (RGS). 

Emission lines from H-like and He-like low-Z ions (from C to Si) and Fe L-shell ions dominate the spec- 
trum. Strong, narrow radiative recombination continua (RRC) for several ions are also present, implying 
that most of the observed soft-X-ray emission arises in low-temperature plasma (fcTe ~ few eV). This 
plasma is photoionized by the inferred nuclear continuum (obscured along our line of sight), as expected 
in the unified model of active galactic nuclei (AGN). We find excess emission (compared with pure re- 
combination) in all resonance lines (Is— >np) up to the photoelectric edge, demonstrating the importance 
of photoexcitation as well. Wc introduce a simple model of a cone of plasma irradiated by the nuclear 
^ . continuum; the line emission we observe along our line of sight perpendicular to the cone is produced 

' through recombination/radiative cascade following photoionization and radiative decay following pho- 

Q> , toexcitation. A remarkably good fit is obtained to the H-like and He-like ionic line series, with inferred 

' radial ionic column densities consistent with recent observations of warm absorbers in Seyfert 1 galax- 

I ies. Previous Chandra imaging revealed a large (extending out to ~ 500 pc) ionization cone containing 

' most of the X-ray flux, implying that the warm absorber in NGC 1068 is a large-scale outflow. To 

I explain the ionic column densities, a broad, flat distribution in the logarithm of the ionization parameter 

. {£_ — Lx/n^r^) is necessary, spanning log^ — 0-3. This suggests either radially-stratifled ionization 

' zones, the existence of a broad density distribution (spanning a few orders of magnitude) at each radius, 

'p^, or some combination of both. 
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that possibility, and was instead c onsistent with recom - 
bination in a photoionized plasma ( Kallman ct al. 1996 ). 
Sako et al. (2000b) pointed out that the same nuclear con- 
tinuum responsible for photoionizing the plasma may also 
photoexcite the ionic resonance transitions, producing line 
emission through radiative decay back to the ground state 
(Krolik & Kriss 1995). Therefore, they asserted that the 
soft-X-ray emission from Markarian 3 could be explained 
in its entirety as reemission in warm plasma both pho- 
toionized and photoexcited by the inferred nuclear contin- 
uum. This model is consistent with reprocessed emission 
from the warm absorbing medium typically observed in 
Seyfert 1 galaxies, but now viewed from a different angle. 

In this paper, we present the first high-resolution X-ray 
spectrum of the prototypical Seyfert 2 galaxy NGC 1068. 
NGC 1068 has a rich observational history, starting with its 
original, thoug h less commonly used, designation as M77 
(Messier 1850). Various independent studies in the optical 
early on confirmed its "nebular" -type emiss ion- line spec- 
trum ( |Fath 1909| ; ^lipher 191% [Hubble 1926|). It appeared 
as the first object in Seyfert's historic list (Seyfert 1942) 



and, more recently, provided t he first convincing evidenc e 
for a unified model of AGN ( Antonucci fc Miller 1985 ). 
Extensive studies at all wavelengths have upheld this 
model: broad optical emission lines in polarized light 



(Antonucci & Miller 1985), the existence of highly coUi- 
mated bipolar radio jets (Wilson & Ulvestad 1982; Wilson 
& Ulvestad 198 3; Pedlar et al. 1 989), biconical O III emis- 



sion regions ( Ebstcin et al. 1989 ), reflected X-ray emission 
of an obscured nucleus (Monier & Halpern 1987; Elvis & 
Lawrence 1988; Koyama et al. 1989), and a soft-X-ray 
component that lies significantly above the refiected con- 
tinuum (MarshaU et al. 1993; Ueno et al. 1994). 

We show below — using a new, simple radiative trans- 
fer code — that the soft-X-ray spectrum of NGC 1068, 
obtained with the XMM-Newton Refiection Grating Spec- 
trometer (RGS), can be explained in its entirety as emis- 
sion from warm, irradiated plasma (as in Markarian 3), 
which reprocesses the nuclear continuum incident on it 
through recombination/radiative cascade following pho- 
toionization and radiative decay following photoexcita- 
tion. 

The format of our paper is as follows. In §^ we begin 
with a discussion of the observations and data analysis. 
Detailed spectral analysis is presented in followed by 
interpretation in terms of possible emission mechanisms in 
§^ In <|^, we discuss a novel method for inferring geomet- 
rical/dynamical properties of the ionization cone from the 
observed relative amount of photoionization versus pho- 
toexcitation. We discuss the remarkably consistent fits we 
obtain using this method in Finally, in we discuss 
the implications of these fits, including the connection be- 
tween the observed large-scale ionization cone and generic 
warm absorbers. 

2. OBSERVATIONS AND DATA REDUCTION 



NGC 1068 was observed with XMM-Newton for 110 ks 
starting at 5:48:05 UT on 2000 July 29. Data were taken 
simultaneously for the EPIC detectors MOS and PN, the 
optical monitor, and the RGS. All data were processed 
with the XMM-Newton Science Analysis Software, SAS 5.2, 
with the corresponding calibration files available for that 
version. We concentrate in this paper mainly on the data 
from the RGS, which was operated in its normal spec- 
troscopy mode. However, in order to confirm the negli- 
gible effect of source size and position on the extracted 
RGS spectrum, we also present a cursory analysis of the 
PN image below. 

The RGS covers the wavelength range of approximately 
6 to 38 A (i; = 0.35-2.5 keV) with a resolution of 0.05 A, 
and a peak effective area of about 140 cm^ at 15 A. The 
SAS filters the observed RGS events in dispersion chan- 
nel versus CCD-pulse-height space to separate the spec- 
tral orders. The background is estimated using events 
from a region spatially offset from the source. The wave- 
lengths assigned to the dispersion channels are based on 
the pointing and geometry of the telescope with an overall 
accuracy of 8 mA (1 -a) across the total wavelength range, 
which is 120 km s~^ at 20 A. Similarly, the estimated 
uncertainty in line broadening is 5 mA (l-cr), which is 
75 km s^i at 20 A. Based on ground calibration, we ex- 
pect the l-fj uncertainty in the effective area to be less 
than 10% above 9 A and at most 20% for shorter wave- 



lengths dden Herder et al. 200l| [den Herder ct al. 2002] ). 

Independent, fiuxed spectra from the RGS 1 (red) and 
RGS 2 (blue) instruments for all 110 ks of data for the 
m — —1 order are presented in Fig. |l[ The spectral dis- 
continuities are due to chip gaps, bad pixels, and the pre- 
vious in-fiight loss of one CCD for RGS 2 (A ~ 20-24 A). 
Overall, the fiuxed spectra agree to <^ 10%. 

3. SPECTRAL ANALYSIS 

3.1. Line Emission 

The soft-X-ray spectrum of NGC 1068 (Fig. |l|) is dom- 
inated by line emission. Emission lines from H-like and 
He-like C, N, O, Ne, Mg, and Si are all clearly detected. 
Numerous Fe L-shell emission lines (Fe XVII to Fe XXIV) 
are present as well. Many higher-order resonance transi- 
tions (Is^np) in H-like and He-like ions labelled Ly/?-e 
and He/?-e, respectively, are prominent, with evidence for 
strong emission from even higher order transitions as well. 
Several unidentified features located in the bottom panel 
of Fig. |l| (e.g., at 27.92, 30.4, 34.0-34.6, and 36.38 A) are 
likely due to L-shell emission from mid-Z ions such as Si 
or S. A bright fluorescent line for Fe K is observed in the 
MOS/PN CCD-resolution data, but the RGS does not ex- 
tend in wavelength down to Fe K. We report marginal 
detection of the fluorescent line of neutral Si in the RGS 
spectrum; there is no other evidence for any other sig- 
nificant fluorescent line emission. We see no significant 
continuum emission in the spectrum. 
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3.2. Line Table 

By performing joint fits to the s eparate m = —1 
and m = — 2 orders within XSPEC (Arnaud 1996), we 
have measured hne fluxes, widths, and shifts for all rel- 
atively bright lines which have been unambiguously iden- 
tified (Table The measured fluxes take into account 
t he Galactic column dens ity of A^^^' = 3.5 x 10^° cm^^ 
(Dickey & Lockman 1990) using the neutral absorption 
model thahs (Wilms, Allen, & McCray 2000) in XSPEC. 
Line shifts are measured with respect to the systemic red- 
shift of NGC 1068 of z = 0.00379 ± 0.00001 or cz = 
1137 ± 3 km s-i (Huchra, Vogeley, & Geller 1999). 

The brightness of NGC 1068 coupled with the duration 
of the observation, implies that all the listed errors (repre- 
senting ±1-ct) of the unblended lines are dominated by in- 
strumental uncertainty (see §^ and not by counting statis- 
tics. For example, there are approximately 6050 and 1880 
total counts in the O VII forbidden line (RGS 1 alone) and 
the Ne X Lya line (RGS 1 and RGS 2), respectively, im- 
plying counting errors of 1.3% and 2.3%; these errors are 
much less than the effective area uncertainty of 10%. 

3.3. Line Shifts and Widths 

We find that most emission lines in the RGS spec- 
trum are shifted and appear broader than expected for 
monochromatic lines (Table p. 

The line shifts are due to velocity shifts at the source 
and not, for example, due to an absolute pointing error. 
The position of the nominal centroid of the X-ray emission 
in the FN image with a CCD-determined energy cut con- 
sistent with the RGS waveband (Fig. ^) has been used to 
derive the RGS spectrum. Gaussian fits to the PN image 
using several different CCD-resolution energy cuts yield 
centroids all within 0".5 of each other, proving that pos- 
sible spatial stratification of different emission zones is in- 
significant. The overall positional uncertainty between the 
PN and RGS instruments is less than the absolute pointing 
uncertainty of the telescope of 4", implying a robust upper 
limit on the line shift uncertainty of 8 mA (120 km s~^ at 
20 A). 

Significant velocity blueshifts from roughly to 
600 km s~^ have been measured from the RGS spectrum 
(Table |l|). Line blueshifts of each ionic line series are fairly 
consistent. Longer wavelength lines appear to have higher 
blueshifts. Overall, the observed blueshifts are comparable 
to blueshifts of optical/UV emission lines (Grimes, Kriss, 
& Espey 1999; Kraemer & Crenshaw 2000), though, in 
contrast to these observations, we find no evidence for any 
significant redshifts. 

The widths we measure are similarly due to intrinsic 
velocity distributions at the source. A Chandra image of 
NGC 1068 (Young, Wilson, & Shopbell 2001) shows evi- 
dence for weak extended emission on scales of tens of arc- 
seconds (Fig. H). The effect of spatial broadening in the 
RGS spectrum for this source, characterized by a strongly 
peaked central component and weak extended emission, 



is negligible (A. Rasmussen, private comm.). Therefore, 
any excess broadening of particular lines can only be due 
to velocity distributions at the source. Measured veloc- 
ity widths are also presented in Table |l|. All widths lie in 
the range 17°'^*' = 300-700 km s~^, which is consistent with 
the linewidths associated with the narrow-line-emitting re- 
gions observed in the UV (Grimes, Kriss, & Espey 1999). 

3.4. Radiative Recombination Continua 

The spectrum also includes very distinctive radiative 
recombination continua (RRC) for H-like and He-like C, 
N, and O, which are produced when electrons recom- 
bine directly to the ground state of these highly-ionized 
species. RRC are broad features for hot, coUisionally- 
ionized plasma, but are narrow, prominent features for 
cooler photoionized plasma. The narrow width of these 
RRC provide a direct measure of the recombining el ectron 
temperature ( Licdahl fc Paerels 1996 ; Liedahl 1999| ). 

For those RRC which are clearly detected (C V, C VI, 
N VI, N VII, O VII, and O VIII), we were able to de- 
termine accurate temperatures from the RRC profile itself 
(Table ||). The unblended RRC (on the short wavelength 
side) of C V and O VII provide the best temperature deter- 
minations. For the other RRC, the dominant source of un- 
certainty comes from blending with other nearby lines. All 
temperatures are consistent with the conservative range of 
fcTe ~ few eV. These temperatures imply velocity broad- 
ening of all features by ^ 10 km s^^. This broadening is 
much lower than the observed broadening of ~ 400 km s~^ , 
implying the importance of bulk and/or turbulent cloud 
velocities. 

Higher temperatures of kT^ = 20-30 eV were inferred 
from the broad, polarized, electron-scattered optical lines 
(Miller, Goodrich, & Mathews 1991). But these tempera- 
tures are really only upper limits, since nonthermal broad- 
ening due to bulk and/or turbulent cloud velocities, as 
observed in the X-ray, could also be present. 

Weak features at the positions of the RRC from Ne IX 
and Ne X are also detected, though severe blending makes 
a temperature determination difficult in these cases. Their 
inferred temperatures are not necessarily discrepant with 
the range given for the other ions, though significantly 
hotter temperatures cannot be excluded. 

3.5. Higher- order Series Lines 

In addition to the prominent n = 2 ^ 1 transitions 
in the H-like and He- like ions, the unique combination of 
high resolution and large effective area characteristic of the 
RGS allows for the detection and flux measurement of sev- 
eral discrete higher-order resonance transitions (np— >ls), 
as well as blended higher-order transitions leading up to 
the RRC edge. For the H-like ions, we list the observed ra- 
tios of these transitions to the Lya line in Table ^, where 
we also give the ratios expected for recombination in a 
warm, photoionized plasma and for hotter plasma in col- 
lisional equilibrium. For the He-like ions, we list the ob- 
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served ratios of these transitions to the forbidden Une in 
Table together with the ratios expected for recombi- 
nation. Here, we do not provide the ratios expected for a 
hotter plasma in coUisional equilibrium, because the strong 
dependence of the relative strength of the forbidden line 
on temperature would make that comparison meaningless. 

3.6. He-like Triplet Ratios 

For He- like N and O, we present measurements of the 
standard plasma diagnostic ratio s R = f/i and G = 
(/ + ^)/'' ( pabriel fc Jordan 1969 ) in Table [g, where r, i, 
and / are the resonance, intercombination, and forbidden 
line fluxes. Also listed are the ratios expected for recom- 
bination in a warm, photoionized plasma and for hotter 
plasma in coUisional equilibrium. 

4. INTERPRETATION OF THE OBSERVED SPECTRAL 
FEATURES 

Unambiguous evidence for the predominance of emis- 
sion in a photoionized plasma comes primarily from the 
observed bright RRC, which indicate recombination at a 
temperature around a few eV. In addition, as we explain 
in detail below, the ratios of the intercombination and for- 
bidden lines to the RRC are very c lose to the values ex- 
pected for pure recombination (see § |6.2| and bottom panel 
of Fig. |l^). The R ratio for N VI (Table |) also appears 
to favor pure recombination, though the R ratio for O VII 
is somewhat ambiguous (see §5.2). 

The anomalous He-like G-ratios (Table ||), however, 
are grossly inconsistent with pure recombination, sug- 
gesting an additional component and/or mechanism. An 
additional, weak collisionally- ionized plasma component, 
whose He-like triplets generically exhibit a higher reso- 
nance line intensity compared to the intercombination and 
forbidden lines, could explain these effects. However, such 
a component would not only enhance the resonance line, 
but a/nines in the triplet, thereby overpredicting the inter- 
combination and forbidden line fluxes relative to the RRC. 
This interpretation is marginally ruled out by the triplet 
to RRC ratio in the RGS data. 

Similarly ambiguous G ratios were previously observed 
in the Chandra HETGS spectrum of the Seyfert 2, Markar- 
ian 3. Sako et al. (2000b) argued against the presence 
of an additional coUisionally-ionized plasma component 
due to the observed weakness of the Fe L shell emission, 
which was consistent with recombination in a photoionized 
plasma (Kallman et al. 1996). Instead, Sako et al. (2000b) 
suggested that the enhanced resonance lines were due to 
photoexcitation (see Fig. ^ for a schematic of the rele- 
vant atomic processes). This argument is sensitive to the 
Fe abundance; however, Sako et al. (2000b) noted that a 
rather unlikely order-of-magnitude underabundance of Fe 
would be necessary to be consistent with the presence of 
additional hot coUisionally-ionized plasma. 

For NGC 1068, we propose a cleaner test based solely 
on each particular ionic line series. The key to ro- 



bustly distinguishing between photoexcitation and addi- 
tional coUisionally-ionized plasma components lies in the 
higher-order series line strengths. We illustrate this new 
method in Fig. ^. Starting with a pure recombination 
spectrum, we add in photoexcitation (self-consistently, as 
presented below in §|5|) or a separate coUisionally-ionized 
plasma component. The addition of photoexcitation en- 
hances all resonance transitions, whereas the addition of 
coUisionally-ionized plasma primarily enhances only the 
triplet lines (with more weight on the resonance line of 
the triplet). 

The ratios of higher-order-series lines to Lya for C VI, 
N VII, and O VIII shown in Table || are significantly higher 
than those expected for coUisional ionization equilibrium 
(CIE), though they are not entirely consistent with pure 
recombination (REC) either. Photoexcitation boosts all 
lines in this series including Lya, but its effect is depen- 
dent on the column density and velocity distribution, so it 
must be modeled carefully. 

However, the forbidden line in He-like ions is unaffected 
by photoexcitation, allowing for a measure of the signifi- 
cance of photoexcitation of higher-order-series transitions 
with respect to the forbidden line, which is fed through re- 
combination after photoionization. Comparison of the ra- 
tios of these higher-ordcr-series lines to the forbidden line 
(Table ^ with similar ratios for plasma in coUisional equi- 
librium is not as useful due to the significant temperature 
sensitivity of /. However, ratios for plasma in coUisional 
equilibrium are generally comparable to or lower than ra- 
tios for pure recombination (Fig. ^). Therefore, the excess 
strength of the higher-order lines, especially the He7 and 
HeJ, can only be explained by photoexcitation. 

5. NOVEL DIAGNOSTICS FROM PHOTOIONIZATION AND 
PHOTOEXCITATION RATES 

We now show that the relative amount of photoioniza- 
tion versus photoexcitation inferred from the soft-X-ray 
spectrum of NGC 1068 provides detailed information about 
the geometry of absorbing/reemitting material around the 
central, power-law-continuum source. 

In the unified model of active galaxies, the observ- 
able properties of a particular AGN depend primarily on 
orientation. The soft-X-ray warm absorber observed in 
Seyfert 1 galaxies, therefore, should show itself in reemis- 
sion for the roughly_perpendicular line of sight to Seyfert 2 
galaxies (see Fig. For the Seyfert 1 view, the rela- 
tive amount of edge absorption versus line absorption de- 
pends greatly on the column density through the material. 
Thus, the Seyfert 2 spectrum should similarly be sensi- 
tive to the radial column density, since each photoioniza- 
tion is balanced by recombination/radiative cascade, and 
each photoexcitation is similarly balanced by radiative de- 
cay. Therefore, the shape of the reemitted Seyfert 2 ionic 
line series can be used to infer the radial (or Seyfert- 1- 
observed) ionic column densities through the warm ab- 
sorber. 
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For low radial ionic column density, N[^^, all lines and 
photoelectric edges remain unsaturated. For a H-like or 
He-like ion, the summed oscillator strength for transitions 
to discrete levels is roughly comparable to the integrated 
oscillator strength for transitions to the continuum (e.g., 
Bethe & Salpeter 1977), implying roughly similar amounts 
of radiative decay following photocxcitation and recombi- 
nation/radiative cascade following photoionization (mod- 
ulated by the shape of the continuum at the relevant pho- 
tocxcitation / photoionization energies) . 

For higher N[^^, the low- lying, high-oscillator-strength 
transitions are the first to saturate, leading to a "tilt- 
ing" of the photocxcitation contribution to the Seyfert 2 
spectrum to lower-oscillator-strength, higher-order-series 
transitions. This can simply be thought of as the stan- 
dard "curve of growth," except now observed indirectly 
through reemission. Due to this gradual line saturation 
and the fact that the edge always saturates last, it is ev- 
ident that an increase in N[^^ implies an increase in the 
overall relative importance of photoionization compared to 
photocxcitation. 

At very high N^^^, all discrete transitions are fully satu- 
rated and the reemitted spectrum is that of pure recombi- 
nation following photoionization. In the following discus- 
sions, these qualitative arguments are made more explicit. 

5.1. Radial Ionic Column Density 

The dependence of our model ionic spectra on radial 
ionic column density, A^j^o^, is illustrated in Fig. 0. Here 
we show the effect of varying N[^^ for He-like O VII. 

The Seyfert 1 and Seyfert 2 views are shown in the top 
three panels on the left and the right of Fig. |^, respec- 
tively, for different radial ionic column densities. The to- 
tal edge absorptions (photoionizations) and line absorp- 
tions (photoexcitations) in the Seyfert 1 view yield the 
expected amounts of recombination/radiative cascade fol- 
lowing photoionization and radiative decay following pho- 
tocxcitation in the Seyfert 2 view. Note the relative 
strength of the lines to the RRC and the effect of sat- 
uration on the high-oscillator-strength, lower-order-series 
transitions. In particular, note the strong variation in the 
relative strength of the resonance line in the He- like triplet 
as a function of optical depth, illustrating that the triplet, 
by itself, is a poor diagnostic for determining whether or 
not hotter, collision-driven plasma (which would show a 
strong resonance line) is present (see Fig. 

5.2. Radial Velocity Width 

The radial velocity width specified by crj,'"^ is a model pa- 
rameter. There is little a priori reason to believe that crj,'*'^ 
should be identical to the observed velocity width of the 
lines, a°^^, viewed perpendicular to the cone (see Fig. ||). 
For this reason, we keep a'if^'^ a free parameter. The effect 
of varying al^"^ is illustrated in Fig. ^. Increasing tr™'^ by 
a certain factor leads to a roughly similar increase (de- 

^http: / /xmm. astro. columbia.edu/research. html 



pending on the line saturation or position on the "curve of 
growth") in all photoexcited lines. This is in contrast to 
variations in N^^^f, which primarily lead to a "tilting" of 
the photocxcitation contribution, since, one-by-one, each 
line, in descending order of oscillator strength, is eventu- 
ally saturated. The velocity width we invoke here may be 
due to a turbulent velocity in the plasma or simply due to 
different bulk velocities of multiple absorption systems. 



6. FITS USING AN IRRADIATED CONE MODEL 

We have constructed a self-consistent model of an ir- 
radiated cone of plasma for application to the NGC 1068 
soft-X-ray spectrum. This model applies to all significant 
H-like and He-like ions, employing ato mic data obtained 
from the Flexible Atomic Code (FAC) ( |Gu 20021 ). (A de- 
tailed study of the more complicated Fe L-shell transitions 
is left for the future.) To simplify fitting, we have incorpo- 
rated our mod el into XSPEC (Arnaud 1996 ) as the local 
model photo'' ( Kinkhabwala et al. 2002| ^ 

In our model, we allow for individual line absorption, 
but assume that photoelectric absorption is negligible. 
With this assumption, we can simplify the calculations 
by irradiating each ionic column density separately with 
an initially unabsorbed power law (discrete line absorp- 
tion due to all other ions can safely be neglected, since it 
removes flux from only a small portion of the total spec- 
trum). Alternatively, using the method employed by our 
model, we can irradiate all ionic column densities together, 
assuming all ions have the same radial distribution. We 
have checked that both methods produce the same repro- 
cessed spectrum as long as all photoelectric edges have low 
optical depths. (The radial ionic column densities we infer 
from the spectrum of NGC 1068 are indeed consistent with 
the assumption of negligible photoelectric absorption, as 
we argue below.) Furthermore, we assume t hat t he plasma 
is optically thin to reemitted photons (see §7.3). 

Our model consists of three free parameters for each ion, 
the radial velocity width crj^'^'^, the covering factor times nu- 
clear luminosity fLx, and the ionic column density N[o^. 
The covering factor, expressed in terms of the solid an- 
gle subtended by the plasma, is / = ft /An, and the total 
power-law luminosity (where, throughout, we assume rea- 
sonable values for the index of F = —1.7 and energy range 
of 13.6 eV-100 keV) is Lx- In order to convert observed 
flux to lumin osity, we take the distance t o NGC 1068 to 
be 14.4 Mpc ( |Bland-Hawthorn et al. 1997 ). 

We assume the following values for the radial velocity 
200 km s^^ and the overall normaliza- 

^43 



width of cr^^d 

tion of fLx = lO"*"^ e rgs s~^ are the same for all the ions 
(discussed further in 6.3). The individual ionic column 
densities, N^^^, are left free for each ion. 

The RRC temperatures are directly observed (Table ^). 
Also, the velocity widths of all emission lines transverse to 
the cone are directly observed (Table |l]). For simplicity. 
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we assume a velocity width of a°^^ ~ 400 km s""'^, which 
works well for all of the observed lines. 

A more complete discussion of the atomic calculations 
and astrophysical assumptions underlying our model is 
presented in Kinkhabwala ct al. (2002). 

6.1. Ftt to H-like CVI 

In Fig. H, we show our fit to the line series correspond- 
ing to H-like C VI. In the bottom panel, we show the best 
fit possible assuming only recombination. Note the gross 
underprediction of all features relative to the RRC. In the 
top panel, we self-consistently add in the contribution from 
photoexcitation assuming the model described above and 
a column density of C VI of Nl^^ = 9 x lO" cm'^. The 
addition of photoexcitation completely explains the excess 
line strengths. 

6.2. Fit to He-like O VII 

Similarly, in Fig. we show our fit to the line series 
corresponding to He-like O VII. In the bottom panel, we 
show the best fit possible assuming only recombination. 
Note the marked underprediction of the resonance line and 
all higher order features relative to the RRC. We note that 
the ratio R = f /i given in Table || is greater than expected 
for pure recombination. This may arise from an additional 
contribution to the forbidden line due to inner-shell pho- 
toionization in O VI (e.g., Kinkhabwala et al. 2002). In 
the top panel, we self-consistently add in the contribution 
from photoexcitation assuming the model described above 
and a column density of O VII of N[^^ = 1.1 x 10^* cm^^ 
Again, all the highcr-order-series lines are completely ex- 
plained by the addition of photoexcitation. 

6.3. Final Fit to All H-like and He-like Ion Series 

Our fit to all H-like and He-like ions in the spectrum 
as a whole (neglecting Fe L shell emission) is presented in 
Fig. |ri|. By "fit," we simply mean a possible set of param- 
eters that can reproduce the spectrum in detail. Though 
photo is a local model in XSPEC, the interdependent, mul- 
ticomponent nature of fitting with this model does not 
all for well-defined error bars for our fitted parameters. 
The best-fit global values for ctJ;^^ and fLx of 200 km s^^ 
and 10*^ ergs s~^, respectively, were found by trial and 
error. Attempts were made to fit the spectrum using 

crj;ad ^ ;^QQ g-l ^rad ^ 4QQ g-1 ^j^j^ g^^^^^j 

malization fLx and individual column densities left free. 
However, the details of the higher-order transitions for any 
given ion were less satisfactorily fit. An estimate for the 
model parameter uncertainty, therefore, gives roughly a 
factor of two for all inferred model parameters: fLx, o'l^'^, 
and N[^^. Ratios of individual column densities, however, 
are much better determined (ranging from a few percent 
to a few tens of percent). Of course, parameters such as 
a°^^ and the RRC temperatures are directly observed with 
error bars quoted in Tables |^ and ||. 



Assuming al""^ = 200 km s~^ and fLx = 10^^ ergs s~\ 
as well as the observed broadening of a°^^ = 400 km s~^, 
we provide all other best-fit parameters in Table ^. 
For each H-like and He-like ion series, all features 
(from the n = 2 ^ 1 transitions up to the RRC) 
are reproduced very well (Fig. ^l]) using these pa- 
rameters. The column densities we infer are indeed 
consistent with our assumption of negligible photoelec- 
tric absorption and, additionally, are consistent (within 
a factor of a few) with the column densities derived 
from several recent high-resolution soft-X-ray spectra of 
Seyfert 1 galaxies (Kaastra et al. 2000 ; Kaspi et al. 200C ; 



Branduardi-Raymont et al. 2001; 3ako et al. 2001). The 
inferred radial velocity distribution of cr™'^ = 200 km s~^ 
is also consistent with the range (crjj^'^ — 100--600 km s~^) 
observed in soft-X-ray spectra of these same Seyfert 1 
galaxies. Further inspection of Table ^ reveals that the 
H-like and He-like ions of each element have approxi- 
mately equal column densities; this important observation 
is discussed further in §7.6, Assuming Solar abundances 
(Grevesse, Noels, & Sauval 1996; Allende Prieto, Lam- 
bert, & Asplund 2001) and a fractional ionic abundance 
of O VIII of /i = 0.5, we o btai n a column density in H of 
iVn ^ 4 X 10^1 cm-2 (see §p]). 



We also give measurements of the emission measure 
of the recombining plasma (EM — J nldV) for each ion 
in Table ^. For all EM values, we assume that the re- 
combining ion fractional abundance is near maximal at 
fi+i = 0.5. For C, N, and O, these emission measures 
are based on the measured RRC temperatures and fluxes. 
For Ne, Mg, and Si, the EM is based on the arbitrarily 
chosen temperature of kT^ = 4 eV. We note that the tem- 
peratures for these ions are not likely to be lower than 
kTc = 4 eV . bu t could easily be an order of magnitude 
higher (see § 7.4 ) , implying the EMs we list are really only 
lower limits for these ions. 

7. DISCUSSION 

We have shown that the soft-X-ray spectrum of the 
prototypical Seyfert 2 galaxy NGC 1068 is dominated by 
emission from warm plasma irradiated by the nuclear con- 
tinuum. In addition, we have shown that the relative 
amounts of photoionization versus photoexcitation for all 
ionic transition series provide detailed, self-consistent in- 
formation about the geometry and dynamics of the emis- 
sion regions, yielding radial column densities, radial veloc- 
ity widths, and the product of covering factor and nuclear 
luminosity. The values we obtain for these geometrical 
properties (in particular, the radial ionic column densi- 
ties) are remarkably similar to absorption measurements 
in Seyfert 1 galaxies. We have also shown that there 
is no measurable additional component due to hot, col- 
lisional plasma (e.g., from shocks, a circumnuclear star- 
burst, or a hot confining medium). We can put a robust 
upper limit to an additional hot, coUisional plasma com- 
ponent in NGC 1068 at an order-of-magnitude less than 
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the line luminosities measured in the soft-X-ray spectrum, 
or EM = / n^dV < 3 x 10^^ cm~'^ for temperatures in 
the range kT — 400-900 eV. In the following, we explore 
further relevant astrophysical considerations. 



7.1. 



Radial Filling Factor of the Warm Absorbing 
Plasma 



The Chandra image of NGC 1068 shown in Fig. 12 



veals that a large fraction of the soft-X-ray emission from 
the ionization cone comes from a region close to the nu- 
cleus, as inferred from radio observations. The minimum 
radius of the ionization cone is clearly < 100 pc, but could 
be much smaller. It is possible that a large fraction of 
the cone is obscured from our view; however, the observed 
photoexcitation of low-n, high-oscillator-strength transi- 
tions implies that we must be directly viewing the bulk 
of the ionization cone, since these high-oscillator-strength 
transitions saturate very quickly (see Fig. |^). This sets a 
lower limit on the radius of the bulk of the ionization cone 
which must be greater than the radius of the obscured 
circumnuclear region. 

Using upper limits to the covering factor (/ < 0.1) and 
minimum radius (rmin 100 pc), we can derive an upper 
limit to the radial filling factor of the plasma g, which, by 
definition, is < 1. We can estimate the column density for 
a given ion as follows: 



ion 



Azf^ 



1.2 



g{r)^dr 



1.6 X 



/^min.pc L0.5J L 10^"^ ergs s ^] [ 1 crgs cm s ^] 1 (1) 

where Az is the elemental abundance, fi is the frac- 
tional ionic abundance, rmin.pc is the minimum radius 
of the cone measured in parsec, g{r) is the radial filling 
as a function of radius, and 5 is a single number rep- 
resenting the average radial filling in the cone. Assum- 
ing the relevant values for C V of Az — 3.63 x 10~^ 
(Grevesse et al. 1996), £, = 1 erg cm s^^, and observed 
column density of Nl^j^'^ = 8 x 10^^ cm^^, we obtain 

ion " 

g = 1.4 X 10~'^/rmin,pc- Similarly, assuming the relevant 
values for O VIII ofAz = 4.89 x 10""* (AUcnde Prieto et al. 
2001), ^ = 10 ergs cm s~^, and observed column density of 
N^^^ = 1 X IQi^ cm-2, we obtain g ^ 1.3 x lO-V^min.pc- 
Using robust upper limits for the covering factor and the 
minimum radius from the Chandra image (Fig. of 
/ < 0.1 and Tmin.pc ^ 100 pc, we obtain upper limits 
oi g < 1.4 X 10^^ (plasma at ^ = 1 erg cm s~^) and 
5 < 1.3 X 10^^ (plasma at ^ = 10 ergs cm s~^). The ac- 
tual filling factor, therefore, is significantly less than unity, 
suggesting that the X-ray plasma fills only a small portion 
of the ionization cone. 

The number densities for plasma with ionization param- 
eter ^ = 1 erg cm s^^ at 500 pc (maximum extent of cone 
from Fig. |l^) down to 1 pc (arbitrarily-assumed minimum 
radius), taking Lx = 10^^ ergs s~-^, are nc = 4 x 10^- 
1 X 10^ cm~'^. These densities are a factor of ten less for 
plasma with ^ = 10 ergs cm s~^. 



7.2. Lack of Electron- scattered Continuum and Nn 

We observe no significant electron-scattered continuum. 
This allows us to obtain an upper limit to the electron 
column density in the context of our model as follows. 
From our measured value for the covering factor times nu- 
clear luminosity of fLx = 10*'^ ergs s~^ and the assumed 
distance to NGC 1068 of 14.4 Mpc, we can write the col- 
umn density in electrons as a function of the normaliza- 
tion A of the reflected continuum FrcR = AE~^ '^ (-Frcfl 
has units of photons cm^^ s^^ keV~^ with E in keV) 
as Nc = 1.0 X IQ-^^A cm~^. The observed upper hmit 
to the electron-reflected continuum of ^ < 8 x lO^"' (us- 
ing the line-free region between 20.0-20.6 A) gives Nc < 



X 10 cm 



We have already f ound a reasonable lower 
-2 



limit to A^H of 4 X 10 cm in §6.3. Therefore, assum- 



ing H and He are fully stripped, we constrain the neutral 
column density to be 4x 10^^ cm~^ < < 7x 10^^ cm~^. 



7.3. Optical Depth in Transverse Direction 

The optical depth in a given line for an approximately 
spherical cloud determines the average number of scatters 
a line photon will experience before it is either destroyed 
or escapes. For significant column density transverse to 
the cone (i.e., along the Seyfert 2 line of sight in Fig. ^), 
the higher-order resonance transitions (np— >ls with n > 2) 
and photons created by recombination directly to ground 
will be reprocessed and degraded to n = 2 ^ 1 transitions 
( |Bowcn 1935| ). 

Since we do not know the geometry of the scattering re- 
gions (e.g., filled cone or isolated clouds), we employ order- 
of-magnitude estimates. If A^j^'^"^ is the ionic column den- 
sity and CTavc is the line-profile- averaged cross section for a 
particular transition (taking into account velocity broad- 
ening ay) in a roughly spherical cloud, then the optical 



thickness is r 



cloud 



(Tavc and the average number of 



scatters is A^ w r. Applying these considerations to H-like 
and He-like series, we find that conversions of Ly/3 to Lya 
and He/3 to Hea provide the most sensitive diagnostics. 
The branching ratios are roughly 0.9 to produce the same 
photon and therefore 0.1 to cascade differently. 

Our final fit to the spectrum of NGC 1068 shown in 
Fig. ^ assumes that all reemission occurs in optically-thin 
plasma, i.e., all photons created in the plasma escape with 
no further reprocessing. Therefore, due to the goodness of 
the fit, we assume at most a 10% reduction (roughly the 
limit of detectability) in the Ly/3 and He/? lines, we plot 
the upper limit to the ionic optical depth versus ay for 
several different ions in Fig. 03. Robust lower and upper 
limits to ay of 10 km s^^ and 500 km s^^ are set by thermal 
broadening [kT ~ 3 eV) and the observed line broadening, 
respectively. He-like ions provide the best upper limits to 
ionic column densities. Notice that the upper limits (over 
the entire range in ay) to the ionic column densities of 
C V, N VI, and O VII are consistent-with-to-lower-than 
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the inferred column densities along the cone given in Ta- 
ble 1^ (where we have assumed a radial velocity distribution 
of al^'^ — 200 km s~^). This may suggest that the radial 
column density is larger by roughly a factor of a few than 
the transverse column density, which would be consistent 
with a narrow cone geometry. 

7.4. Electron Temperatures 

The electron temperatures we measure directly from 
the RRC for recombination onto specific ions can be com- 
pared with results from self-consistent simulations of pho- 
toionized plasma, where photoelectric heating is balanced 
by line cooling in the plasma. Using the photoioniza- 



tion code XSTAR ([Kallman fc Krolik 1995| ), we have de- 
termined the expected temperature, kT, versus ioniza- 
tion parameter, ^ = Lx/rier'^, for a spherical, optically- 
thin plasma with central, ionizing, power-law continuum 
(r = -1.7, S,nin = 13.6 eV, E^,^ = 100 keV). In Fig. |l|, 
we present the expected range in temperature and ioniza- 
tion parameter for all observed H-like and He-like ions. 
The XSTAR-predicted temperatures agree well with the 
data. 

7.5. Thermal Instability 

For a generic photoionized plasma, thermal instabilities 
may develop. In Fig. |l^, we plot the thermal equilibrium 
curve for the same optically thin plasma (assuming So- 
lar abundances) along which heating equals cooling. The 
plasma is thermally unstable along segments of the equilib- 
rium curve which have negative slope, since upward devia- 
tions from the curve in this region heat the plasma until it 
reaches the higher stable branch and downward deviations 
cool the plasma until it reaches the lower stable branch. 
Comparing both Figs. ^ and |l^, we see that several ions 
present in the spectrum have maximal ionic abundances 
in this thermally unstable region. This putative thermal 
instability, therefore, appears problematic for ions in the 
observed range of ionization parameter, ^. 

We note that thermal instability can be removed if ad- 
ditional heating/cooling mechanisms are at work in the 
plasma in addition to those associated with photoelectric 
heating and line cooling. Thermal instability is also sen- 
sitive to metal abundance (Hess, Kahn, & Paerels 1997), 
with low metal abundances effectively removing the insta- 
bility. The existence of thermal instability in astrophysical 
photoionized plasmas is presently an open question. 

7.6. Ionization Parameter Distribution 



In Fig. 16, we present the fractional ionic abundances 
relative to H and assuming Solar abundances for the el- 
ements observed in the RGS spectrum of NGC 1068. No 
single ionization parameter is capable of explaining the 
column densities obtained in Table ^, which shows sim- 
ilar column densities in H-likc and He-like ions for each 
elemental species. The data instead require a broad, flat 
distribution in ionization parameter of the plasma over 



the range log^ ^ 0-3. This distribution in ^ = Lx/ricr"^ 
could be due to spatial stratification of io nization zones in 
a single-density cone QSako et al. 2000a| ) or a broad den- 
sity distribution at each radius. 

Kraemer and Crenshaw (2000b) find that the extended 
optical/UV absorber is powered by photoionization due to 
an inferred hidden source of continuum, as we similarly 
argue for the soft- X-ray emission. These findings, cou- 
pled with the remarkable overlap of the soft-X-ray image 
of the cone with optical observations of O III emission 
regions (Young et al. 2001), favor an intrinsic density dis- 
tribution (over at least three orders of magnitude) at each 
radius. For such a model, a perfectly flat distribution in 
log ^ would give a distribution in densities at each radius of 
/(no) oc n~^. An intrinsic density distribution in AGN ion- 
ization cones has been proposed in the past to explain the 
optical and UV spectra of generic Seyfert 2 galaxies (e.g., 
Komossa & Schulz 1997). In X-ray emitting plasmas, such 
a density distribution may naturally arise in outflows due 
to thermal instabilities (Krolik & Kriss 2001), which may 
counter the usual perception of such instabilities o n th e 
presence of specific ions in the un stab le region (see §7.?;). 

The upper limits obtained in §7J_ for the radial filling 
factor do not allow for distinguishing between spatially- 
stratified ionization zones or an intrinsic density distribu- 
tion. Robust discrimination between these two models will 
require future observations and analysis. 

7.7. Relative Elemental Abundances 

Given the ionization parameter distribution, the in- 
ferred ionic column densities can provide estimates of rel- 
ative elemental abundances. However, even absent knowl- 
edge of the exact distribution, some tentative conclusions 
can be drawn due to the overlap in ionization parameter 
of specific charge states of different ions. 

The necessity of a broad, rather flat range of ioniza- 
tion parameter allows for comparison of overlapping charge 
states of different elements, such as H-like C and He-like O 
(see Fig. ^6|). Overall, we find that all of the well- measured 
ionic column densities yield relative elemental abundances 
which are consistent with Solar (Grevesse et al. 1996; Al- 
lende Prieto et al. 2001), excluding N, which is overabun- 
dant by approximately a factor of three. The column den- 
sities of Mg and Si also seem high in the context of a flat 
ionization parameter distribution. But, due to low counts 
and blending, their column densities are not well deter- 
mined, making difficult any robust conclusions concerning 
their relative abundance. 

Furthermore, we note that the abnormally high 
Fe/0 abundance invoked for previous ASCA spectra of 
NGC 1068 by Netzer & Turner (1997) is not observed. In 
fact, the O VII and O VIII emission lines are the strongest 
in the spectrum, especially when compared to Fe L-shell 
transitions, which are considerably weaker (but not nec- 
essarily inconsistent with Solar- abundance photoionized 
plasma). Fe L-shell line emission in generic photoionized 



9 



plasma, as well as in NGC 1068, is currently under inves- 
tigation. 

Wc would like to thank Julian Krolik and the anony- 
mous referee for helpful comments. AK acknowledges use- 
ful discussions with all members of the Columbia XMM- 
Newton RGS group. This work is based on observations ob- 
tained with XMM-Newton, an ESA science mission with in- 
struments and contributions directly funded by ESA Mem- 
ber States and the USA (NASA). The Columbia University 
team is supported by NASA. AK acknowledges additional 



support from an NSF Graduate Research Fellowship and 
NASA GSRP fellowship. MS and MFG were partially sup- 
ported by NASA through Chandra Postdoctoral Fellow- 
ship Award Numbers PFOl-20016 and PFOl-10014, respec- 
tively, issued by the Chandra X-ray Observatory Center, 
which is operated by the Smithsonian Astrophysical Ob- 
servatory for and behalf of NASA Lmdcr contract NAS8- 
39073. SRON is supported by the Netherlands Founda- 
tion for Scientific Research (NWO). Work at LLNL was 
performed under the auspices of the U. S. Department of 
Energy, Contract No. W-7405-Eng-48. 



References 

AUende Prieto, C, Lambert, D. L., Asplund, M. 2001, ApJ, 556, 1, L63 
Antonucci, R. R. J. 1993, ARA&A, 31, 473 
Antonucci, R. R. J., & Miller, J. S. 1985, ApJ, 297, 621 

Arnaud, K. A. 1996, in ASP Conf. Ser. 101, Astronomical Data Analysis Software Systems V, ed. G. H. Jacoby & J. 
Barnes (San Francisco: ASP), 17 

Bethe, H. A. & Salpeter, E. E. 1977, Quantum Mechanics of One- and Two-Electron Atoms (New York: Plenum 
Publishing Corporation) 

Bland-Hawthorn, J., Gallimore, J. F., Tacconi, L. J., Brinks, E., Baum, S. A., Antonucci, R. R. J., Cecil, G. N. 1997, 
Ap&SS, 248, Issue 1/2, 9 

Bowen, I. S. 1935, ApJ, 81, 1 

Branduardi-Raymont, G., Sako, M., Kahn, S. M., Brinkman, A. C, Kaastra, J. S., & Page, M. 2001, A&A, 365, L140 
den Herder, J. W., et al. 2001, A&A, 365, L7 

den Herder, J. W., et al. 2002, to appear in the proceedings of the symposium 'New Visions of the X-ray Universe in the 
XMM-Newton and Chandra Era', 26-30 November 2001, ESTEC, The Netherlands 

Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 

Ebstein, S. M., Carleton, N. P., & Papaliolios, C. 1989, ApJ, 336, 103 

Elvis, M., & Lawrence, A. 1988, ApJ, 331, 161 

Path, E. A. 1909, Lick Observatory Bulletin, 5, 17 

Gabriel, A. H., & Jordan, C. 1969, MNRAS, 145, 241. 

Grevesse, N., Noels, A., & Sauval, A. J. 1996, in ASP Conf. Ser. 99, Cosmic Abundances, ed. S. S. Holt & G. Sonneborn 
(San Francisco: ASP), 117 

Grimes, J. P., Kriss, G. A., & Espey, B. R. 1999, ApJ, 526, 130 
Gu, M. F. 2002, in preparation 

Hess, C. J., Kahn, S. M., & Paerels, F. B. S. 1997, ApJ, 478, 94 
Hubble, E. 1926, ApJ, 54, 369 

Huchra, J. P., Vogeley, M. S., & Geller, M. J. 1999, ApJS, 121, 287 

Kaastra, J. S., Mewe, R., Liedahl, D. A., Komossa, S., & Brinkman, A. C. 2000, A&A, 354, L83 

Kallman, T. R., & Krolik, J. H. 1995, XSTAR = A Spectral Analysis Tool, HEASARC (NASA/GSFC, Greenbelt) 



10 

Kallman, T. R., Liedahl, D., Osterheld, A., Goldstein, W., & Kahn, S. 1996, ApJ, 465, 994 

Kaspi, S., Brandt, W. N., Netzer, H., Sambruna, R., Chartas, G., Garmire, G. R, & Nousek, J. A. 2000, ApJ, 535, L17 
Kinkhabwala, A., Behar, E., Sako, M., Gu, M. F., Kahn, S. M., & Paerels, F. 2002, in preparation 
Komossa, S., & Schulz, H. 1997, A&A, 323, 31 

Koyama, K., Inoue, H., Tanaka, Y., Awaki, H., Takano, S., Ohashi, T., & Matsuoka, M. 1989, PASJ, 41, 731 

Kraemer, S. B., & Crenshaw, D. M. 2000a, ApJ, 532, 256 

Kraemer, S. B., & Crenshaw, D. M. 20G0b, ApJ, 544, 763 

Krolik, J. H., & Kriss, G. A. 1995, ApJ, 447, 512 

Krolik, J. H., McKcc, C. P., & Tarter, C. B. 1981, ApJ, 249, 422 

Liedahl, D. A. 1999, in X-Ray Spectroscopy in Astrophysics, ed. J. van Paradijs & J. A. M. Bleeker (Berlin: Springer), 
189 

Liedahl, D. A., & Paerels, F. 1996, ApJ, 468, L33 

Marshall, F. E., ct al. 1993, ApJ, 405, 168 

Messier, C. 1850, Connaissancc dcs Temps, 1784, 227 

Miller, J. S., & Antomicci, R. R. J. 1983, ApJ, 271, L7 

Miller, J. S., Goodrich, R. W., & Mathews, W. G. 1991, ApJ, 378, 47 

Monier, R., & Halpern, J. P. 1987, ApJ, 315, L17 

Muxlow, T. W. B., Pedlar, A., Holloway, A. J., Gallimore, J. F., & Antonucci, R. R. J. 1996, MNRAS, 278, 854 
Netzer, H., & Turner, T. J. 1997, ApJ, 488, 694 

Pedlar, A., Booler, R. V., Spencer, R. E., & Stewart, O. J. 1983 MNRAS, 202, 647 
Sako, M., Kahn, S. M., Paerels, F., & Liedahl, D. A. 2000a, ApJ, 542, 684 
Sako, M., Kahn, S. M., Paerels, F., & Liedahl, D. A. 2000b, ApJ, 543, L115 
Sako, M., et al. 2001, A&A, 365, L168 
Seyfert, C. K. 1943, ApJ, 97, 28 

Sliphcr, V. M. 1917, Lowell Observatory Bulletin, 3, 59 

Smith, R., & Brickhouse, N. 2000, Astrophysical Plasmas: Codes, Models, and Observations, Proceedings of the con- 
ference held in Mexico City, October 25-29, 1999, ed. J. Arthur, N. Brickhouse, & J. Franco, Revista Mexicana de 
Astronomia y Astrofisica (Serie de Conferencias), 9, 134 

Ueno, S., Mushotzky, R. F., Koyama, K., Iwasawa, K., Awaki, H., & Hayashi, I. 1994, PASJ, 46, L71 

Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 

Wilson, A. S., & Ulvestad, J. S. 1982, ApJ, 263, 576 

Wilson, A. S., & Ulvestad, J. S. 1983, ApJ, 275, 8 

Young, A. J., Wilson, A. S., & Shopbell, P. L. 2001, ApJ, 556, 6 



11 




oo 



CD 



o 



CD 



oo 

CN 



[l_V 3_lu3 suo^oqd ^_ol] xny 



Fig. 1.— Effective-area-corrected, first-order RGS 1 (red) and RGS 2 (blue) spectra of NGC 1068 shifted to its rest frame (z = 0.00379). 
The spectral discontinuities are due to chip gaps in the CCD arrays, bad pixels, and the previous in-flight loss of one CCD for RGS 2 
(A ~ 20—24 A). Line labelling indicates the final state ion. All H-like (a) and He-like (r, i, and f) principal order lines are labelled for each ion 
with the corresponding RRC edges indicated as well. Additionally, resonance transitions (np— »ls) are labelled as (3 through S (short for Ly/3-(5 
and He/3-5). Several Fe L-shell transitions are listed as well. Unlabelled features at longer wavelengths (e.g., A = 27.45,27.92,30.4,31.0, 
34.0—34.6, and 36.38 A) are likely due to L-shell tranisitions in mid-Z elements. Line blueshifts are especially noticeable at longer wavelength. 
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Fig. 2.-- Image of NGC 1068 taken with the PN camera on board XMM-Newton. Only X-rays with CCD-determined energies less 
than 2.5 kcV were taken. We use factor-of-two contours with the lowest at 2.8 X 10^* counts arscc^'^ s^^. The image is centered at the 
peak of emission at a(J2000)=02h42™40=.65, <5(J2000)=-00°00'41".40 (1" = 72 pc). The RGS dispersion axis is indicated (oriented 72.34° 
counterclockwise from the vertical), with A indicating the direction of increasing wavelength. 



60 




-60 I , , , \ , 1 1 I 1 1 iJlj 1 1 : I , , , \ , , 1 

60 40 20 -20 -40 -60 

RA direction [arcsec] 



Fig. 3. — Chandra ACIS-S contour image of NGC 1068 using a frame time of 0.4 s (Young et al. 2001). Only X-rays with CCD-determined 
energies less than 2.5 keV were taken. We use factor-of-ten contours with the lowest at 3.5 x 10"* counts arsec"^ s~^. The image is centered 
on the peak of the X-ray emission at a(J2000)=02''42™40=.71, (5(J2000)=-00°00'47".9 (1" = 72 pc). 
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Fig. 4. — Atomic processes of photoexcitation and photoionization with inverse processes of radiative decay and recombination/radiative 
cascade, respectively. Most radiative decays after photoexcitation will occur directly to the ground state, but other paths (as shown in the 
diagram) are also possible. 
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Fig. 5. — We demonstrate how to distinguish between hot coUisional plasma (e.g., starburst region) and photoexcitation with the He-like 
O VII line series. Starting with pure recombination following photoionization ("PI" — top two panels, fcTe = 4 eV), we self-consistently 
add radiative decay following photoexcitation ("PE"- middle left panel), assuming reasonable ionization cone parameters (as in §B|), or an 
additional hot plasma component in coUisional ionization equilibrium ("CIE" — middle right panel, kT = 150 oV), for which we obtain the 
bottom two panels. Note that both bottom panels have similar triplet ratios, implying that consideration of the triplet alone is insufficient 
to discriminate between the two scenarios. However, the "AGN ALONE" panel has significantly stronger higher-order-series transitions 
(including the RRC) than the "AGN-I-STARBURST" panel, demonstrating the diagnostic importance of these transitions. Photoexcitation 
generically enhances higher-order-series transitions beyond the amount possible in a hot, collisionally-ionized plasma due to the flatness of the 
incident photoionizing spectrum (compEired with a Maxwellian electron distribution) as well as the characteristic "curve of growth" effects 
obtained in irradiated plasmas (see Fig. M). (Normalization in each panel is arbitrary.) 
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Seyfert 1 

A f=^l/47t 




Fig. 6. — Simple cartoon of ionized cone model (not to scale). The central nuclear component, denoting the black-hole/accretion- 
disc/comptonizcd-halo system, is shown as the black spot. In the Seyfert 1 view, the ionized cone (warm absorber) is seen through its 
absorption of the intrinsic nuclear continuum Lx ■ In the Seyfert 2 view, the intrinsic continuum is highly absorbed by intervening material 
("dusty torus" drawn to the right and left of the nucleus), allowing for the reprocessed emission in the ionized cone to be observed. The 
variables comprising the two global parameters a]^'^ and fLx in our model are indicated. The radial ionic column density N^^^, which is 
left free for each ion, is also indicated. The Seyfert- 2-observed velocity distribution perpendicular to the cone is rj'^^ ■ For NGC 1068, the 
intrinsic continuum is likely completely obscured, contributing no flux to the soft-X-ray regime. Also, in NGC 1068, the NE cone is much 
brighter than its counterpart in the SW (hence the asymmetry in the above diagram), therefore, the covering factor / applies to the NE cone 
alone. Generic Seyfert 2 galaxies, however, may show evidence for two equally bright ionization cones. 
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Fig. 7. — Effect of differing radial column densities on the reemitted spectrum for He-like O VII. The top three panels on the left show 
the radial Seyfert 1 view through the outflow and down to the nucleus for radial column densities in O VII of 10^®, 10^^, and 10^^ cm~^. 
(The logarithm of the O VII cross section for photoexcitation and photoionization with separating boundary is presented in the lower left 
panel.) Corresponding panelsjDn the right show the Seyfert 2 view roughly perpendicular to the axis of outflow and from which the nucleus 
is completely obscured (Fig. a). All photons absorbed out of the power law in the Seyfert 1 spectrum are reprocessed and reemitted to 
generate the Seyfert 2 spectrum. Radiative decay following photoexcitation dominates the Seyfert 2 spectrum at low column densities, 
whereas recombination/radiative cascade following photoionization dominates at high column densities (for comparison, pure recombination 
is shown in the lower right panel). For H-like ions, the behavior is similar (omitting the intercombination/for bidden lines). Throughout we 
take a radial gaussian distribution with cr^^'* = 200 km s~^, Seyfert-2-observed velocity distribution with <J°^^ = 400 km s~^, and temperature 
kTc = 3 eV. (Normalization in each panel is arbitrary.) 




Fig. 8. — Effect of different crj^^"^ for He-like O VII spectrum using a radial column density of 10^® cm ^. At larger crj^^'^, photoexcitation 
is enhanced relative to photoionization. Spectra on the right-hand side were convolved with the same transverse velocity distribution 
= 400 km s-i) and assume kTa = 3 eV. 
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Fig. 9. — Final fit parameters for C VI (see Table |^ and Fig. |ll| ) including recombination/radiative cascade following photoionization and 
radiative decay following photoexcitation (top). Recombination alone (bottom) is unable to explain the excess emission in all resonance lines 
np— >ls. (Data are plotted in rest frame of NGC 1068. The model spectrum has been velocity shifted by —400 km s~^.) 
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Fig. 10. — Final fit parameters for O VII (see Table ^ and Fig. including recombination/radiative cascade following photoionization 
and radiative decay following photoexcitation (top). Recombination alone (bottom) is unable to explain the excess emission in all resonance 
lines np— >ls. (Data are plotted in rest frame of NGC 1068. The model spectrum has been velocity shifted by —400 km s~^.) 
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Fig. 11. — Data from RGS 1 (blue) plus the best fit model (red) to NGC 1068, which assumes a normalization of fLx = Iff*^ ergs s~^ 
and radial velocity distribution of crjj^'* = 200 km s^^, with individual ionic column densities given in Table H. We assume a total neutral 
column density to the source of Nn = 3.5 X 10^" cm~^ and an observed velocity broadening of a-°^^ = 400 km s~^. The data have been been 
shifted to the rest frame of NGC 1068 (z = 0.00379). The model wavelengths are not shifted. Line blueshifts, especially noticeable at long 
wavelengths, are clearly present. Fe L shell emission (spanning 9 A to 18 A) has not yet been included in the model. 
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Fig. 12. — Chandra ACIS-S contour image of NGC 1068 using a frame time of 0.4 s and total exposure of 11465.4 s (Young et al. 2001). 
Only X-rays with CCD-determined energies less than 2.5 keV were taken. The image is centered on the peak of the X-ray emission at 
a(J2000)=02'i42™40^71, (5(J2000)=-00°00'47".9 (1" = 72 pc). The position of th e peak of the 5 GHz radio emission (assumed nuclear 
position) at a(J2000)=02''42™40=.715, (5(J2000)=-00°00'47".64 ( |Mux1ow ct al. 1996| ) is within 0".3 of the X-ray peak, which is well within 
the Chandra positional uncertainty. We use factor-of-two contours with the lowest at 4.4 X 10~^ counts arsec"^ s~^. The cone is clearly 
extended over hundreds of parsec. 




Fig. 13. — Upper limits to ionic column densities assuming 10% conversion percentage of Ly/3 photons to Lyo and, similarly, He/3 to 
Hea due to multiple scatterings in a roughly isotropic medium. Depending on the filling factor, the medium can be thought of as either an 
individual cloud or as the entire emission region. (In order to further generalize these results for arbitrary line conversion, each line in the 
plot can be multiplied by In// In 0.9, where 1 — / is the degradation of the specific Sp^ls transition.) The labels in the upper left and lower 
right reflect the order of the lines, from top to bottom, in the plot. Robust upper and lower limits to (T„ of 10 km s~^ and 500 km s~^ are 
from thermal broadening (kT = 3 eV) and the observed emission line broadening, respectively. For H-like C VI and He-like C V|J^ VI, and 
O VII, the upper limits over the entire range are consistent with or significantly less than the fitted radial column densities (Table fel), possibly 
suggesting a narrow cone geometry. 
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Fig. 14. — From a photoionization simulation using XSTAR of an optically-thin, low-density cloud of Solar-abundance plasma surrounding 
an ionizing point source with power-law spectrum (F = 1.7, -Emin = 13.6 eV, -Emax = 100 keV), we have determined the expected relationship 
between temperature kT [eV] and ionization parameter ^ = Lx /(nor'^) (solid curve), where Lx [ergs s~^] is the total luminosity in the 
power-law continuum, rie [cm~^] is the electron density, and r [cm] is the distance from the nucleus to the emission region. Peaks in the 
recombination emissivity forming H-like and He-like ions are indicated (dotted lines). The FWHM for each ion's emissivity distribution is 
A(log^) 1. The x's and o's denote the temperature confidence interval and best fit, respectively, taken directly from Table H. 




Fig. 15. — For the same simulation as in Fig. we have determined the curve of thermal stability using H = ^/47rfcT (Krolik, McKee, & 
Tarter 1981). Regions along the thermal equilibrmm curve that have negative slope are thermally unstable (cooling is greater than heating 
for regions leftward oLthe curve and vice versa for regions rightward). Several ions observed in the spectrum of NGC 1068 have maximal 
abundance (see Fig. h4) in the thermally unstable region. 
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Fig. 16. — Fractional ionic abundances relative to H (assumin^Solar ahundances) for He-like (leftmost curves), H-like (middle), and bare 
(rightmost) ions using the same_XSTAR simulation as in Figs. h4 and ha. The observed similar column densities in the H-like and He-like 
charge states of each ion (Table p|) require a broad distribution in ^. Assuming relative Solar abundances all ions, the observation of column 
densities in C V and C VI of 8—9x10^^ cm~^ and in Si XHI and Si XIV of 2 X 10^^ cm~^ requires this distribution to be rather flat as well, 
spanning the range log^ = 0-3. 
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N VII Lya 


24.781 


-270 ± 100 


650 ± 60 


6.09 ±0.74 


N VI He5 


23.277 


-210 ± 100 


470 zt 60 


0.32 ±0.14 


N VI He7 


23.771 


-490 zt 100 


360 zt 60 


0.77 ±0.12 


N VI r 


28.787 


-350zb80 


410 zt 50^ 


3.84 ±0.38 


N VI i 


29.083 


-310zb80 


410 ± SO'' 


1.22 ±0.15 


N VI f 


29.534 


-430 zt 80 


410 zt 50^ 


8.46 ± 0.85 


C VI Ly6 


26.357 


-440 zt 90 


400 zt 60 


0.57 ±0.10 


C VI Ly7 


26.990 


-440 ± 90 


320 zt 60 


0.95 ±0.11 


C VI Ly/3 


28.466 


-380 zt 80 


450 zt 50 


2.03 ±0.20 


C VI Lya 


33.736 


-360 zt 70 


510 zt 40 


12.29 ± 1.23 


C VHe(5 


32.754 


-560 zt 70 


510 ±50 


1.33 ±0.20 


C VHei 


33.426 




550 zt 50 


1.50 ±0.31 


CVHe/J 


34.973 


-550 zt 70 


360 zt 40 


1.28 ±0.30 



"^Error bars (due to dominant systematic instrumental uncertainty) are 1-cr. 

''Fe L-shell lines and He-like triplet linewidths for N and O were tied together for fitting for each line complex. 



Table 2 



Ion RRC kT,. FliLx [Ur'^ pliolojis cm-- 



C V 


9 c+0.5 


4.30 ± 0.43 


C VI 


4.0^};° 


2.83 ±0.28 


N VI 




2.06 ±0.21 


N VII 


4.0t^;0 


1.14±0.11 


VII 




2.43 ±0.24 


VIII 




1.25 ±0.13 



''Error bars denote approximate 90% confidence interval. 
''Error bars (due to dominant effective- area uncertainty) are 1-a. 



Table 3 



Ion 


Line Ratio 


CIE^ 




Observed'^ 


C VI 


Ly/3/LyQ: 


0.090 


0.144 


0.165 ±0.023 




Ly7/LyQ: 


0.024 


0.052 


0.077 ±0.012 




LyJ/Lya 


0.010 


0.026 


0.046 ± 0.009 


N VII 


Ly/3/Lya 


0.096 


0.141 


0.133 ±0.025 




Ly7/Lya 


0.023 


0.049 


0.071 ±0.014 




LyS/Lya 


0.009 


0.024 


0.020 ±0.010 


VIII 


Ly/3/Lya 


0.100 


0.138 


< 0.247 ±0.035 




Ly7/LyQ; 


0.025 


0.048 






LyS/Lya 


0.010 


0.023 


< 0.063 ±0.013 



^Evaluated at CIE temperature which maximizes ion fraction. 

''Evaluated at kT^ = 4 eV, although the temperature dependence 
is weak. 

■^Error bars are 1-a. 



Table 4 



Ion 


Line Ratio 


REC* 


Observed'^ 


N VI 


He/?// 


0.047 


< 0.12 ±0.06 




Hc7// 


0.016 


0.091 ±0.017 




Re5/f 


0.008 


0.038 ±0.017 


VII 


He/?// 


0.051 


0.138 ±0.019 




Hc7// 


0.017 


0.083 ±0.012 






0.008 


0.065 ±0.010 



^Evaluated at kT^ = 3 eV for N VI and 4 eV for 
O VII. 

''Error bars are 1-a. 
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Table 5 



Ion 




Ratio 




REC'" 


Observed'^ 


N VI 




R = !/i 


5.0-5.5 


6.0-7.0 


6.9 ± 1.1 


VII 






4.0-4.3 


3.4-4.0 


5.5 ±0.8 


N VI 


G 


^if + i)/r 


0.5-1.0 


4.0-5.5 


2.5 ±0.2 


VII 


G 




0.5-1.1 


4.0-5.3 


2.5 ±0.2 



'^Valid over realistic r ange of coUisional kTc for each ion 
([Smith fc Brickhouse 2000| ). 



•^Vahd for kT^ = 0.1-100 eV ( [Kinkhabwala et al. 200^ . 
'^Error bars are 1-a. 



Table 6 



Ion 


7V-f [cm-2] 


kTe [eV]" 


EM [10'^'' cm-3] 


C V 


8E17 


2.5 


6.7 


C VI 


9E17 


4.0 


1.8 


N VI 


6E17 


3.0 


6.5 


N VII 


6E17 


4.0 


1.7 


VII 


1.1E18 


4.0 


1.5 


VIII 


1E18 


4.0 


0.32 


Ne IX 


3E17 


4 


0.42 


NeX 


2.5E17 


4 


0.097 


Mg XI 


2E17 


4 


0.31 


Mg XII 


2E17 


4 


0.092 


Si XIII 


2E17 


4 


0.13 


Si XIV 


2E17 


4 


0.042 



^kTc values in italics are arbitrarily assumed. 

''Assumes Solar abundance and fractional ionic abundance 
/i+i = 0.5. Values in italics correspond to italicized (and arbi- 
trary) kTe values. 



